function y = func_n2_from_p(p, CONSTS)
    eps = CONSTS.eps;
    g = CONSTS.g;
    q2 = func_q2_from_p(p, CONSTS);
    
    y = -(eps./(p*g)).*(p.^2+q2.^2+(g^2/eps)-eps);
end
